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Resumen 


Este trabajo está integrado en tres partes, en la primera se citan brevemente ideas 
generales sobre las sequías, hasta llegar a su definición y tipos. En la segunda parte 
se expone con detalle el procedimiento sugerido para la estimación de las curvas 
de severidad-área-frecuencia (SAF), y en la tercera parte se describe la aplicación 
numérica al valle de la ciudad de San Luis Potosí, México. La secuencia propuesta 
consta de las cinco etapas básicas siguientes: (1) recopilación y selección de registros 
autilizar; (2) interpolación espacial de registros; (3) cálculo de los déficits anuales de 
lluvia; (4) procesamiento probabilístico de los déficits anuales, y (5) construcción de 
las curvas SAF. Las características temporales y espaciales de las sequías son básicas 
para entender su comportamiento y evolución dentro de una cuenca o región 
específica. Usando uno de los índices más simples de caracterización de las sequías 
y una técnica elemental de ponderado espacial, se establece el procedimiento de 
estimación regional de las curvas SAE, el cual utiliza la distribución General de 
Valores Extremos para obtener las predicciones que definen la frecuencia de tales 
curvas. Por último, se formulan las conclusiones, las cuales destacan la utilidad de 
las curvas SAFE, para la comparación de las sequías de una región con otra y para la 
interpretación de los años de sequía histórica. También se manifiesta la sencillez del 
método propuesto y se indica que su primera etapa es la más crítica. 


Palabras clave: precipitación mensual, valores medianos, homogeneidad, inter- 
polación espacial, predicciones, distribución GVE. 


Introducción 


Con el objeto de afrontar de manera racional 
la escasez de agua, es necesario entender 
los procesos que gobiernan su ocurrencia y 
sus causas. En general, la escasez de agua 
tiene dos orígenes: el primero es natural y se 
debe a la aridez y a las sequías; el segundo es 
inducido por las actividades humanas, por ello 
incluye a la desertificación y a sus demandas o 
requerimientos hídricos. Las sequías son uno 
de los desastres naturales más dañinos que 
ocurren casi cada año en alguna región del 
mundo; causan, directa o indirectamente, daños 
económicos, sociales y ambientales. Las sequías 
son inevitables, pero sus impactos pueden ser 
mitigados (Mishra y Singh, 2008; Tallaksen y 
Van Lanen, 2004). 


Las sequías son provocadas por los cuatro 
mecanismos siguientes: (1) sistema océano- 
atmósfera, (2) anomalías de la temperatura 
superficial del océano, (3) relaciones tiempo 
atmosférico-actividad solar, y (4) procesos de 
origen de los huracanes y sus anomalías. Debi- 
do a que los mecanismos anteriores no actúan 
de manera independiente sino combinada, 
resulta sumamente difícil el estudio y la 
predicción de las sequías. De hecho, las sequías 
son uno de los peligros naturales más complejos 
y menos entendidos (Mishra y Singh, 2008; 
Sheffield y Wood, 2011). 

Aparte de la complejidad en sus orígenes, 
las sequías presentan grandes dificultades 
para cuantificar sus daños, ya que afectan 
a muchas personas y a diferentes sectores 
económicos, además de que su extensión 
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geográfica por lo común es extensa y entonces 
se complica de forma notable la estimación 
de sus impactos ambientales y personales. 
Ante tales dificultades, la investigación sobre 
sequías se ha orientado a establecer índices 
y/o procedimientos que intentan definirlas, 
detectarlas y medirlas (Tallaksen y Van Lanen, 
2004; Sheffield y Wood, 2011). 

Existiendo la necesidad de contrastar las 
sequías de una región a otra, así como también 
comparar las diferentes sequías históricas 
ocurridas en la misma región, Mishra y Singh 
(2008) han propuesto un enfoque de análisis 
similar al regional de frecuencia de crecientes, 
el cual permite estudiar a las sequías en 
términos de su duración, severidad, extensión 
espacial y probabilidad de ocurrencia, para 
construir sus curvas de severidad-área- 
frecuencia (curvas SAF). El objetivo de este 
trabajo consiste en exponer con detalle el 
procedimiento elemental formulado para la 
estimación de las curvas SAE, describiendo su 
aplicación numérica en la cuenca del valle de 
la ciudad de San Luis Potosí, México. 


Estimación de las curvas SAF 
Definición y clasificación de las sequías 


Las sequías difieren de otros peligros naturales, 
como las crecientes, los huracanes y los terre- 
motos, en que sus efectos se acumulan durante 
meses e incluso años; además, es difícil preci- 
sar su inicio y terminación. La definición más 
general indica que es un periodo sostenido 
y prolongado de escasez de agua. Como tal 
deficiencia de agua también está asociada con 
su Origen o naturaleza, se identifican cuatro 
tipos de sequías: (1) la meteorológica, definida 
como la falta de precipitación sobre una región 
durante un cierto lapso; (2) la hidrológica, 
relacionada con un periodo con disponibilidad 
de escurrimiento insuficiente para satisfacer 
las demandas; (3) la agrícola, ligada por lo 
general a una deficiencia de humedad del 
suelo y la consecuente pérdida de los cultivos, 
y (4) la socio-económica, asociada con la falla del 


sistema de recursos hidráulicos que abastece 
las demandas, y entonces origina problemas 
sociales y daños en las actividades económicas. 


Etapas del procedimiento 


La estimación de las curvas SAF para analizar 
sequías meteorológicas implica procesar los 
registros de precipitación mensual con base 
en las siguientes cinco etapas para su esta- 
blecimiento: (1) recopilación y selección de 
registros a utilizar; (2) interpolación espacial de 
registros; (3) cálculo de los déficits anuales de 
lluvia; (4) procesamiento probabilístico de los 
déficits anuales, y (5) construcción de las cur- 
vas SAF. En la figura 1 se muestra el diagrama 
de flujo de las etapas citadas, y se detallan 
los análisis o procedimientos específicos que 
incluyen. 


Recopilación y selección de registros 
pluviométricos 


Esta primera etapa requiere del conocimiento 
geográfico de la región o cuenca para la cual 
se estimarán las curvas SAF. Lógicamente, 
entre más información pluviométrica esté 
disponible, las curvas SAF serán más con- 
fiables; lo mismo ocurre si la zona o área es 
plana y uniforme en condiciones físicas. 
Por el contrario, la escasez de registros o la 
disponibilidad de amplitud reducida con- 
ducen a curvas SAF inexactas. 

La primera parte de esta etapa implica la 
recopilación de toda la información pluvio- 
métrica disponible. En seguida, se toman en 
cuenta sus periodos de registro y se visuali- 
zan los problemas asociados con la deducción 
de todos los datos faltantes, sean éstos men- 
suales o anuales. Lo anterior, orientado a 
definir el mayor lapso de registro común 
en todas las estaciones. Posteriormente, se 
aborda la selección de los registros que serán 
utilizados, para lo cual se toma en cuenta su 
distribución dentro de la región o área, y 
su representatividad, según las condiciones 
geográficas de ésta, para buscar disponer de 
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Etapa del procedimiento: 


Tarea o técnica específica: 


11 Obtención de toda la información pluviométrica disponible. 


12 Contraste de lapsos de los registros disponibles 


(selección exhaustiva de pluviómetros y del periodo común). 


I. Recopilación y selección 
de registros pluviométricos 
por utilizar. 


13 Definición de datos anuales faltantes. 

14 Definición de datos mensuales faltantes. 

I5 Estimación de datos mensuales y anuales faltantes. 

L6 Selección del periodo común. 

1.7 Distribución de lluvia mensual. 

18 Verificación de la homogeneidad estadística de los registros. 


19 Selección y reemplazo de estaciones pluviométricas 


(selección restringida del periodo común y de pluviómetros). 


Il. Interpolación espacial 
de registros. 


11.1 
111.2 
TII.3 


TIT. Cálculo de los déficits 
anuales de la lluvia. 


1v.1 
Iv.2 
IV.3 


IV. Procesamiento 
probabilístico de los 
déficits anuales. 


V. Construcción de las 
curvas SAF. 


A | II.1 Construcción de la malla de cuadrados. 


11.2 Ponderación espacial de registros mensuales. 


Selección del índice de sequía por utilizar. 
Selección de las duraciones de sequía por procesar. 
Aplicación del índice de sequía seleccionado. 


Selección de los periodos de retorno por procesar. 
Selección del modelo probabilístico por utilizar. 
Selección del método de ajuste o de los resultados. 


V.1 Mapas de resultados por duración y periodo de retorno. 
V.2 Proceso de selección de cuadros contiguos máximos. 
V.3. Dibujo de las curvas de severidad-área-frecuencia (SAF). 


Figura 1. Diagrama de flujo de las etapas y los procedimientos que las integran en el método sugerido para la 


construcción de las curvas SAF. 


una red lo más amplia y uniforme posible. 
Despué se realiza la estimación de todos los 
datos faltantes de cada registro seleccionado, 
mensuales y anuales en el periodo común. 
Por último, se lleva a cabo una verificación 
de requerimientos estadísticos de los regis- 
tros seleccionados, con base en sus registros 
anuales. Por lógica, esta etapa concluye con 
una nueva selección y remplazo, cuando algu- 
na delas estaciones, previamente seleccionada, 
no resultó homogénea. 


Interpolación espacial de registros 


Aunque en la etapa anterior se hubiera logrado 
establecer una red de estaciones pluviométricas 
bastante uniforme dentro de la región de 
estudio, sus áreas de influencia de acuerdo 
con los polígonos de Thiessen probablemente 
serían muy pocas para poder tomar en cuenta 
la distribución espacial de las sequías en 
la construcción de las curvas SAF. Por ello, 
Mishra y Singh (2008) recomiendan formar 
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una cuadrícula sobre la región o cuenca, que la 
divida en al menos 25 áreas uniformes, y luego 
crear en el centro de cada cuadro un registro 
hipotético con base en los reales disponibles. El 
razonamiento de soporte para la interpolación 
espacial es el hecho de que los registros de 
precipitación son similares entre más cercanos 
están y difieren más conforme están más 
alejados. Entonces, la técnica básica y elemental 
de interpolación es la ponderación de cada 
registro real con el inverso de la distancia al 
centro de la malla (Mishra y Singh, 2008). 

Ahora cada precipitación mensual hipo- 
tética del periodo común en el centro de 
cada cuadrado de la malla es igual a la suma 
ponderada de la misma precipitación mensual 
en cada registro real. Se acepta que los pesos o 
factores de ponderación (a,) son inversamente 
proporcionales a las distancias (d,) entre cada 
sitio del registro hipotético (Y) y cada uno de 
los m registros reales (X,), esto es (Campos- 
Aranda, 2008): 


m 
Misa =Ya, o con 


3 


a, =1.00 (1) 


1 


—d(1/d, +1/d,+-+1/d,) Q) 


Az 


El índice ¡ es mensual y varía de 1 a 12; 
en cambio, el índice ¡ es anual y por lo tanto 
fluctúa de 1 al número de años (11) del periodo 
común procesado. 


Cálculo de los déficits anuales de lluvia 


Un índice de sequía es una variable que indica 
el efecto de la escasez de agua y que permite 
definir diversas características de cada sequía, 
por ejemplo, duración y severidad. Por su- 
puesto que existen varios índices para la 
sequía meteorológica que han sido sugeridos y 
todos ellos emplean la precipitación sola o en 


combinación con otros elementos climáticos, 
como la temperatura. Uno de los enfoques más 
simples ha sido expuesto por Campos-Aranda 
(2012) y consiste en obtener la suma de las 
deficiencias de la precipitación mensual con 
respecto a su mediana, para obtener el déficit 
anual. La duración de las sequías es tomada 
en cuenta al definir los periodos más lluviosos 
por analizar: (1) un mes, (2) los tres meses 
consecutivos, (3) los seis meses consecutivos y 
(4) todo el año. 

Para el cálculo del déficit anual en cada 
centro de cuadro de la malla, se tienen que 
estimar los valores medianos, contra los que se 
compara la lluvia hipotética que se obtuvo en la 
interpolación espacial, esto se realiza también 
con base en las ecuaciones (1) y (2). 


Procesamiento probabilístico de los déficits 
anuales 


Al tomar en cuenta que las predicciones 
asociadas con los déficits anuales obtenidos en 
cada centro de la cuadrícula, para duraciones 
de 1, 3, 6 y 12 meses, que son necesarias para 
construir las curvas SAE, tendrán periodos de 
retorno relativamente bajos, de 5, 10, 25, 50, 75 
y 100 años, que se pueden obtener mediante el 
ajuste de cualquiera de las distribuciones de 
probabilidad actualmente recomendadas para 
procesar valores máximos, como son la General 
de Valores Extremos, la Log-Pearson tipo Il y 
la Logística Generalizada, pues es probable que 
ellas conduzcan a predicciones similares. 


Integración de las curvas SAF 


Para realizar esta etapa se deben elaborar 24 
mapas de la región o cuenca, a fin de indicar 
en cada centro de su cuadrícula una predicción; 
entonces son cuatro para cada uno de los 
seis periodos de retorno procesados, siendo 
cada uno de ellos el correspondiente a una 
de las duraciones empleadas. Con base en 
cada mapa se calculan los valores máximos 
asociados con áreas crecientes, desde el cuadro 
del valor máximo hasta el promedio de todos 
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los cuadros. Se sigue un procedimiento que 
intenta reproducir el avance de las sequías 
en el mundo real. Se identifica el cuadro con 
la predicción máxima, después se revisan los 
cuadros contiguos a sus lados y se selecciona 
el de valor máximo, se obtiene su promedio y 
se continúan revisando ahora todos los cua- 
dros adyacentes, para hacer la suma y promedio 
de tres cuadros. Se sigue así sucesivamente 
hasta obtener el promedio de todas las pre- 
dicciones indicadas en el mapa procesado, 
el cual equivale al 100% de área (Sheffield y 
Wood, 2011). En la figura 2 se esquematiza por 
pasos la selección de las magnitudes máxi- 
mas, y la definición de sus valores contiguos 
o aledaños. 

Con los resultados del procesamiento de 
cada mapa se construye una curva SAF en un 
gráfico natural, en cuyas abscisas se indican 
las áreas del 0 al 100%, o bien de O al tamaño 
en km? de la región o cuenca procesada. 
En primer punto dibujado corresponde al 


porcentaje o área en km” que tiene un cuadro 
de la malla construida; los puntos posteriores 
son estos valores acumulados para dos, tres y 
así sucesivamente hasta el número de cuadros 
de la malla o valores con que se define cada 
curva SAF. Lógicamente, en las ordenadas se 
dibujan los valores de las predicciones de los 
déficits, obtenidos como máximos por áreas 
contiguas. Se forman cuatro gráficas de curvas 
SAF, una para cada duración analizada, y cada 
gráfica tiene seis curvas, una para cada periodo 
de retorno procesado. 


Utilidad de las curvas SAF 


Como ya se indicó, las curvas SAF caracterizan 
a las sequías meteorológicas de una zona O 
región, y por lo tanto la comparación entre 
ellas permitirá apreciar sus diferencias, 
principalmente en intensidad, es decir, en 
las magnitudes dibujadas en el eje de las 


ordenadas. Las curvas SAF también permiten 


Paso 1) selección de valor máximo. 


Paso 2) máximo y valores contiguos. 


FB 


a 


A yA A 


Figura 2. Definición por pasos de valores máximos por áreas contiguas (Tallaksen y Van Lanen, 2004). 


Tecnología y 


= Ciencias del Agua, vol. IV, núm. 3, julio-agosto de 2013 


Ciencias del Agua, vol. IV, núm. 3, julio-agosto de 2013 


10)) 


Campos-Aranda, Caracterización de sequías meteorológicas mediante curvas de severidad-área-frecuencia, en el valle de San Luis Potosí,... 


Bojouda] 


el contraste de un año específico, para observar 
su frecuencia y evolución espacial en el 
contexto de tal región. 


Aplicación numérica en el valle de San 
Luis Potosí 


Información pluviométrica disponible 


Dentro del valle donde se localiza la ciudad de 
San Luis Potosí (SLP); existen ocho estaciones 


pluviométricas con una distribución bastante 
irregular, como se aprecia en la figura 3. Aunque 
la estación Zaragoza, ubicada hacia el sureste 
de la ciudad, está fuera del parteaguas, debido 
a su cercanía es considerada representativa de 
tal zona y por ello se tomó en cuenta en los 
análisis posteriores. 

Las características generales de estas nueve 
estaciones pluviométricas y de otra ubicada al 
norte y fuera del valle se presentan en el cuadro 
1. Toda la información de lluvias mensuales 


>  Parteaguas del valle de SLP. 
To >: Cauces de arroyos o ríos. 
vw Carreteras estatales y vialidades. 
x cg del cuadrado de la malla. 
a,b,c Estación pluviométrica no utilizada. 
10) Estación pluviométrica utilizada. 
+ + 
101” 15' 101* 00' 


+22" 15 


+ 222 00' 


1 Estación pluviométrica Los Filtros. 

2 Estación pluviométrica El Peaje. 

3 Estación pluviométrica Mexquitic. 

4 Estación pluviométrica El Terrero. 

5 Estación pluviométrica Xoconoxtle. 
Nota: carta topográfica de apoyo F14-4. 


+ 2230 


100" 45' 


70 


100” 45" 


Figura 3. Cuenca del valle de la ciudad de San Luis Potosí y malla de cuadrados para su discretización. 


Campos-Aranda, Caracterización de sequías meteorológicas mediante curvas de severidad-área-frecuencia, en el valle de San Luis Potosí,... 


disponible en las estaciones pluviométricas 
del cuadro 1 fue obtenida en las oficinas de 
la Dirección Estatal San Luis Potosí de la 
Comisión Nacional del Agua (Conagua). 

En la columna final del cuadro 1 se citan los 
años que están faltantes en cada registro y que 
por lo tanto se tendrán que estimar mediante 
regresión con algún registro pluviométrico 
cercano que los contenga. Aparte existen los 
años incompletos, en los cuales faltan de uno a 
cinco meses, mismos que serán estimados con 
base en la moda y la precipitación confiable 
mensuales, siempre y cuando éstos no 
pertenezcan en su totalidad al periodo lluvioso 
de mayo-octubre, pues en tal caso el año fue 
considerado faltante. 

El único registro que no presenta años in- 
completos en su periodo total disponible es el 
del observatorio meteorológico de la Facultad 
de Ingeniería de la Universidad Autónoma 
de San Luis Potosí (UASLP), el cual se puede 
obtener hasta el año 2000 en Moreno-Castro et 
al. (2004), y de 2001 a 2009 en Campos-Aranda 
(2012). Las lluvias mensuales del año 2010 son 
las siguientes: 27.9, 84.7, 0.0, 14.7, 29.9, 31.1, 
53.4, 32.6, 0.0, 0.9, 4.0 y 0.0. Por el contrario, 
uno de los registros más incompletos es el 
del observatorio meteorológico de San Luis 
Potosí, ubicado en el Parque Tangamanga Il, 


antes el aeropuerto, con ocho años faltantes en 
el periodo 1971-1978. En los cuadros 2 y 3 se 
exponen los registros de la estación El Peaje, 
completa, y de El Terrero, la más incompleta. 
Debido a la cercanía de la estación Zaragoza 
con Xoconoxtle, la primera será procesada para 
contraste de las estimaciones en la segunda. 


Establecimiento del periodo común de datos 


El periodo común para el estudio fue definido 
por el inicio de operación de la estación 
Xoconoxtle, la cual comenzó su registro en 
octubre de 1980; por lo anterior, abarcará de 
1981 a 2010 con sólo treinta años. La estación 
El Terrero, que comenzó en 1983, requerirá 
de la estimación de sus dos años iniciales. 
Aunque la amplitud del periodo común puede 
parecer corta, se considera adecuada debido a 
que el periodo de retorno de cien años será el 
más grande procesado para definir el término 
de frecuencia en las curvas SAF. Conviene 
recordar (Kite, 1977) que es frecuente en los 
análisis probabilísticos considerar aceptable la 
extrapolación del registro de tres a cinco veces, 
para obtener la mayor predicción. Todas las 
estimaciones que se describen a continuación, 
relativas a meses y años faltantes, así como 
los diversos análisis que se realizarán con 


Cuadro 1. Características generales y del periodo común 1981-2010 de las estaciones pluviométricas del valle 
de San Luis Potosí y cercanas. 


Coordenadas geográficas Periodo de Años faltantes (núm. de 
Nombre Clave registro en incompletos) entre 
Long. Latitud Altitud* años 1981-2010 
1. Facultad de Ingeniería - 101* 01” 22" 08' 1893 1922-2010 - 
2. Los Filtros (DGE) 24069 101* 01” DOO9 1877 1949-2010 - 3 
3. Observ. SLP (SMN) 24111 100* 59* 22911 1877 1966-2010 - 3 
4. El Peaje 24024 101* 06' DORDOS 2 100 1963-2010 - E 
5. Xoconoxtle 24179 100* 42 237 05 2 035 1980-2010 (7) 7 
6. Mexquitic 24042 101* 07' DOSMLOS 2 062 1942-2010 1993, 1994 (5) E 
7. Soledad Diez Gtz. 24081 100% 55' 22" 11 1 882 1945-2010 1984, 1989, 1990, 1998 (6) E 
8. El Terrero 24197 I0OH5SÍ 229 24 1890 1983-2010 1981, 1982, 1997, 1998 (11) a 
9. Zaragoza 24106 100* 44” 229 03 2 000 1961-2010 1984 (1) E 
10. Los Pilares 24038 (IO DOEDIS 1778 1963-2008 2009, 2010 S 
* En metros sobre el nivel del mar (msnm). 3 
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Cuadro 2. Precipitación mensual (mm) en el periodo común en la estación pluviométrica El Peaje del valle de 
San Luis Potosí. 


Año Ene. | Feb. | Mar. | Abr. | Mayo | Jun. Jul. Ago. Sep. | Oct. | Nov. | Dic. | Anual 
1981 34.9 | 18.0 25 | 416 76.4 114.0 60.0 25.5 63.0 57.4 10 10.0 504.3 
1982 0.0 1.0 19 27.0 65.0 1.0 SEL 30.0 TUS 59.8 | 18.9 15.8 291.4 
1983 24.5 0.0 0.0 0.0 51.4 41.3 61.5 23.3 38.8 13.0 5.5 0.0 259.3 
1984 57.0 | 16.5 0.0 0.0 30.0 CSul | hs So LOES A ESSi0) 0.0 9.6 507.7 
1985 10 0.0 7.5 28.0 75.8 156.5 99.0 76.7 4.2 18.2 7.7 10.7 | 485.3 
1986 0.0 3.0 0.0 | 90.8 29.8 75.8 89.5 23.0 54.9 74 | 270 2.0 403.2 
1987 13.2 | 10.2 0.0 | 34.0 30.0 62.2 98.7 93 62.3 0.0 259 0.0 322.6 
1988 175 0.0 A 50 225) 42.3 70 | 125 48.7 2.0 0.0 AS 399.4 
1989 2.0 15 0.0 10 37.0 80.4 3.0 100.5 24.0 0.0 55 28.0 282.9 
1990 oO 13.8 UE 3.0 67.0 20.5 321 48.1 724 | 101.8 10 0.0 5727 
1991 6.0 6.0 0.0 0.0 53.8 39.5 | 219.1 19.4 86.2 305 | 155 0.0 476.0 
1992 119.9 | 26.0 9.0 15.9 | 104.4 26.5 SAD 13.3 870 56.7 | 24.7 3 477.1 
1993 17.4 5.0 0.0 7.9 21.5 123.7 | 42.4 12.5 134.8 6.0 33.5 0.0 404.7 
1994 17.0 0.0 0.0 | 44.0 38.7 97.0 40.8 26.3 41.3 40.9 2.0 0.0 348.0 
1995 7.0 0.0 15.0 0.0 63.4 125 24.5 | 138.8 16.9 8.5 10.7 27.4 324.7 
1996 0.5 0.0 0.0 0.0 20.2 46.5 3502 764 | 1000 | 325 0.0 0.0 391.3, 
1997 26.5 | 18.5 43.6 | 67.2 58.0 525 | 118.0 15.0 12 31.5 7.0 0.0 439.0 
1998 0.0 5.0 0.8 0.0 2.8 46.4 12.0 27/03 AO 9S:S 0.0 0.0 339.3 
1999 0.0 0.0 59 5.7 27.7 225 | 131.9 25.8 95.4 9.4 3.7 71 334.7 
2000 0.7 0.0 51 0.0 3907 | 1095 187 38.0 29 21.6 YE 28.8 281.9 
2001 13 3.0 150 | 21,5 18.7 64.0 40.7 53.0 74.9 13.2 2.0 14.5 321.8 
2002 43.0 | 29.8 0.0 10.1 13.4 90.3 91.3 69.5 ME) | 532 | Sell 0.0 556.7 
2003 8.0 0.0 0.0 11.7 5.3 318 | 121.1 29.0 31.5 80.0 0.0 0.0 319.0 
2004 42.4 3.0 88.4 6.5 235 16700 SO S7.0 98.3 46.3 0.0 OS 702.7 
2005 2.3 422 | 11.50 2.8 3.8 46.0 59.4 26.7 59.5 18.7 0.0 0.9 273.8 
2006 13 0.0 0.4 18.3 54.0 15.0 78.2 | 1343 | 147.1 | 56.2 | 26.9 12 556.4 
2007 12.0 | 145 16.9 | 18.7 58.8 120.6 80.3 85.6 82.9 62 | 23.7 0.0 520.2 
2008 0.0 1.8 0.0 | 32.6 18.6 NS | 107 | 2622 80.7 Ja 0.0 0.0 581.1 
2009 0.0 0.0 3.4 3.0 45.8 512 | 1195 38.7 | 175.7 | 69.9 6.0 20.6 533.6 
2010 32.3 | 126.3 0.0 12.0 18.7 629 | 257 | 520 54.1 0.5 0.0 0.0 563.9 

Valor 

ntaleña 0.0 0.0 0.0 0.0 2.8 1.0 3.0 9.5 1.2 0.0 0.0 0.0 259.3 

Valor 

máximo 119.9 | 126.3 88.4 90.8 104.4 167: Z1O 1 725272 175.7 | 101.8 | 38.1 28.8 702.7 
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los registros pluviométricos mensuales com- 
pletados para llegar a construir las curvas SAF 
estarán basados exclusivamente en el periodo 
común 1981-2010. 


Selección por ubicación de las estaciones 
pluviométricas 


Como ya se indicó, la ubicación de las ocho 
estaciones pluviométricas disponibles den- 


tro del valle de San Luis Potosí no es ni 
remotamente ideal. En realidad, la mitad de 
ellas están localizadas al suroeste del valle y 
tres de éstas dentro de la zona urbana; además, 
tales estaciones son las de registro completo en 
el periodo común: Facultad de Ingeniería, Los 
Filtros, Observatorio SLP y El Peaje. De tales 
estaciones, El Peaje, por ubicarse en la zona 
cerril, será utilizada, así como alguna de las 
restantes. Por su proximidad con las estaciones 
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Cuadro 3. Precipitación mensual (mm) en el periodo común en la estación pluviométrica El Terrero del valle de 
San Luis Potosí. 


Año Ene. | Feb. | Mar. | Abr. | Mayo | Jun. Jul. Ago. Sep. Oct. | Nov. | Dic. | Anual 
1981 (1Q2| (6D | (661) | (128) | (35.6) | (49.4) | (73.0) | (34.7) | (51.0) | (16.8) | (8.3) | (13.1) | (313.1) 
1982 (7.9) | (2.0) | (.0) | (8.3) | (23.1) | (32.1) | (47.4) | (22.6) | (33.2) | (10.9) | (5.4) | (8.5) | (203.4) 
1983 (9.8) | (2.5) | (Q.5) | (10.3) | (28.8) | (39.8) | (58.8) | (28.0) | (41.1) | (13.6) | (6.8) | (10.2) | (252.2) 
1984 (0.0) | (0.0) | (0.0) 0.0 25.6 60.0 83.6 22.0 48.0 5.0 0.0 10 | 272 


1985 1.0 0.0 1.0 17.0 84.0 86.0 149.8 14.0 16.0 8.0 2.0 10.0 388.8 
1986 0.0 1.0 0.0 25.0 25.0 739 43.0 8.0 60.0 20 270 4.0 29D) 
1987 8.0 11.0 0.0 15.0 21.0 72.0 72.0 19.0 58.0 0.0 2.0 0.0 278.0 


1988 (13.7) | (8.5) | (8.5) | (14.4) | (40.0) | (55.5) | (81.9) | (39.0) | (57.3) | (18.9) | (9.3) | (14.6) | (351.6) 
1989 (6.3) (1.6) | (1.6) (6.6) | (18.4) | (25.5) | (87.7) | (17.9) | (26.4) | (8.7) | (43) | (6.7) | (161.7) 
1990 (7.6) LS) | (19) (8.0) | (22.3) | (30.9) | (45.6) | (21.7) | (31.8) | (10.5) | (5.3) | (8.2) | (195.7) 
1991 0.0 13.0 (0.0) 0.0 45.0 74.0 124.0 6.0 67.0 28.0 0.0 29.0 386.0 
1992 106.0 3.0 10.0 10.0 46.0 27.0 30.0 8.0 11.0 45.0 | 20.0 4.0 320.0 
1993 60.0 0.0 (0.0) 91.5 38.0 83.0 31.0 18.0 134.0 3.0 | (0.8) | (0.0) 459.3 
1994 16.0 0.0 (0.0) 41.0 19.0 163.0 29.0 35.0 23.0 57.0 2.0 1.0 386.0 
1995 2.0 (0.0) 0.0 3.0 48.0 4.0 46.0 | 216.0 35.0 1.0 11.0 20.0 386.0 
1996 (0.0) (0.0) | (0.0) (145) 65.0 3.0 340 | 114.0 116.0 | 21.0 3.0 0.0 SO7iO 
1997 (11.8) | (3.0) | (3.0) | (12.4) | (34.5) | (47.8) | (70.5) | (33.6) | (49.3) | (16.3) | (8.1) | (12.5) | (302.8) 
1998 (6.4) (1.6) | (1.6) (6.8) | (19.0) | (26.3) | (38.8) | (18.5) | (27.1) | (9.0) | (4.5) | (6.9) | (166.5) 


1999 0.0 0.0 5.0 4.0 2.0 2.0 73.0 38.0 72.0 11.0 3.0 2.0 212.0 
2000 0.0 0.0 0.0 0.0 6.0 71.0 7.0 10.0 0.0 17.0 6.0 8.0 125.0 
2001 0.0 0.0 0.0 22.0 12.0 46.0 13.3 67.0 62.0 2.0 0.0 0.0 224.3 
2002 (0.0) (0.0) 0.0 0.0 14.0 42.0 101.0 | (2.2) (LS) | 7840 | 150 (0.0) 258.5 
2003 8.0 3.0 0.0 4.0 16.0 53.0 61.0 6.0 33.0 37.0 4.0 86.5 311.5 
2004 13.0 0.0 26.6 1.0 58.0 104.0 26.0 (212) (LS) | (22) | (0) 0.0 245.0 
2005 4.0 29.0 6.0 (1.5) 3.0 2.0 60.0 25.0 48.0 0.0 1.0 0.5 180.0 


2006 10.0 0.0 8.0 4.5 63.5 8.0 70.0 56.0 95.0 19.0 | 20.0 14.0 368.0 
2007 14.0 16.0 8.0 (1.5) 62.0 (2.0) 95.0 54.0 50.0 15.0 | 26.0 0.0 343.5 


2008 AO O AA AA AA E AAA E 
2009 00 | 10 | 90 | 270 | 680 | 130 | 505 | 180 | 320 | 190 | 380 | 60.6 | 336.1 
2010 50 lo lso la less ass ima ls lao o loo 00 ase 
Nela 00 | 00 | 00 | 00 | 20 2.0 zo | 22 | oo | o0 | 00 | 00 | 1250 
minimo 
Valor 


E 106.0 | 29.0 26.5 91.5 84.0 163.0 | 149.8 | 216.0 | 134.0 73.0 | 38.0 86.5 459.3 
máximo 


Nota: entre paréntesis, los valores estimados. 


de la zona urbana, Soledad Diez Gutiérrez no además se escogerá una estación de la zona 
fue seleccionada, pero será completada para urbana. 

usarse como estación auxiliar, al estimar los 

valores faltantes de las estaciones Mexquitic y Estimación de valores mensuales faltantes 

El Terrero. En resumen, El Peaje, Mexquitic, El 

Terrero y Xoconoxtle, son seleccionadas debido Se insiste en que sólo se procesan, por meses, 
a su ubicación próxima a la periferia del valle, todos los datos disponibles en el periodo de 
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1981-2010 para estimar los valores mensuales 
de la moda y de la precipitación confiable, definida 
esta última como la lluvia con probabilidad 
de excedencia del 75%. Tales estimaciones 
se realizaron ajustando a cada serie mensual 
disponible del dato faltante la distribución 
Gamma Mixta (Campos-Aranda, 2005). Por 
tratarse de un estudio de sequías, los valores 
buscados deben ser bastante confiables, así se 
aceptó la siguiente regla de selección: cuando 
existe la moda se adopta ésta, siempre y cuando 
no exceda a la precipitación confiable, de ser el 
caso, se adopta tal valor. Lo anterior se ilustra 
en el cuadro 4, con negritas para la estación 
Mexquitic, que es de las más incompletas. 


Distribución mensual de la lluvia 


El fin primordial de la estimación de los valo- 
res mensuales faltantes, a través de su moda o 
de la precipitación confiable, es completar los 
años incompletos y así disponer del mayor 
número de años para obtener magnitudes 
confiables del promedio mensual y de sus 
porcentajes respecto del promedio anual. Estas 
cantidades se tienen en el cuadro 5. Con base en 
tales porcentajes, se distribuirá mensualmente 
cada valor anual estimado por regresión. 


Estimación de valores anuales faltantes 
Cuando en un registro faltó uno o más años 


completos, o bien existen años que presentan 
más de cuatro meses sin datos, y éstos coin- 


ciden en el periodo lluvioso, se estimó su valor 
a través de regresión lineal con la estación 
pluviométrica más cercana y éste se reparte 
según los porcentajes definidos en el cuadro 
5. En el cuadro 1 se observa que existen cuatro 
estaciones pluviométricas con años faltantes: 
Soledad Diez Gutiérrez, Mexquitic, El Terrero 
y Zaragoza. 

El diagrama de dispersión de la lluvia anual 
entre la estación Soledad D.G. y Observatorio 
SLP, que es la más cercana, permitió detectar 
dos datos dispersos: los años 1992 y 1998. 
Al eliminar tales valores, se obtuvo la recta 
de regresión con 24 parejas de datos, cuyo 
coeficiente de correlación lineal fue de 0.805, 
con un error estándar de la estimación de 
54.3 mm (Campos-Aranda, 2010). Los cuatro 
valores anuales estimados se muestran entre 
paréntesis en el cuadro 6. Desafortunadamente, 
los diagramas de dispersión entre Soledad D.G. 
y las estaciones incompletas Mexquitic y El 
Terrero no mostraron correspondencia lineal. 

La estación El Terrero es la más incompleta, 
pues le faltan los dos años iniciales del periodo 
común, así como 1997 y 1998. Las lluvias 
mensuales de 1983 se tomaron como valores 
anormales altos, pues su suma fue de 970.7 
mm. Por el contrario, los valores mensuales de 
los años 1988, 1989 y 1990 fueron considerados 
anómalos bajos, ya que sus sumas anuales 
fueron de 48.7, 16.0 y 22.8 mm, respectivamente, 
entonces estos años también están faltantes. 
En resumen, le faltan ocho años en el periodo 
común. Para completar la estación El Terrero, 


Cuadro 4. Resultados en milímetros de la estimación de la moda y la precipitación confiable en la estación 


pluviométrica Mexquitic con base en la distribución Gamma Mixta. 


Concepto Ene. | Feb. | Mar. | Abr. | Mayo | Jun. | Jul. | Ago. | Sep. | Oct. | Nov. | Dic. 
Número de datos 26 29 28 29 27 28 28 29 28 30 28 27 
Número de valores nulos 11 11 16 13 1 dl 0 2 0 2 10 9 
Media aritmética 14.3 1L5 6.8 137 | 47 | 722 | 791 | 537 | 71.7 | 32.9 6.9 8.1 
Moda NE* NE 4.8 NE YE 23.4 9 19.4 29.6 1.4 2.8 NE 
Mediana muestral 3.4 1.6 0.0 1.2 40.1 | 805 | 748 | 38.0 | 644 | 26.6 4.9 2.0 
Mediana poblacional 319) 2.4 0.0 2.4 ss 00 | Ed2 [222 | 50 | 224 3.8 Sl 
Precipitación confiable 0.0 0.0 0.0 0.0 14.7 27.8 27.0 | 199 32.4 8.6 0.0 0.0 


* NE: no existe. 
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Cuadro 5. Distribución de la precipitación media mensual (mm) en las estaciones pluviométricas del valle de San Luis Potosí. 


oa Ene. | Feb. | Mar. | Abr. | Mayo | Jun. | Jul. | Ago. | Sep. | Oct. | Nov. | Dic. | Anual 
y valor 

Fac. de Ing. 

Lluvia media 141 | 10.0 3.4 19.9 | 38.7 | 649 | 585 | 488 | 70.9 | 28.1 8.9 Y MS YSES, 
% del anual 3.8 2.7 0.9 5.3 10.4 | 174 | 15.7 | 131 | 190 7.5 2.4 2.1 | 100.3 
Los Filtros 

Lluvia media 15.5 9.6 5.1 20.3 | 413 | 61.1 | 789 | 494 | 66.2 | 31.0 8.5 74 | 394.3 
% del anual 39 2.4 ES, Ball 105 | 153 | 200 | 15 11685 YE) 220) 19 | 100.0 
Observ. SLP 

Lluvia media 15.4 95 515 222 | 418 | 63.7 | 646 | 45.3 | 622 | 30.9 7.8 CN 5 70:/ 
% del anual 41 25 15 5.9 11.1 170 | 172 | 121 | 16.6 8.2 2.1 18 | 100.1 
Soledad D. G. 

Lluvia media 149 | 10.1 3.4 223 | 343 | 56.5 | 56.0 | 39.8 | 559 | 27.2 8.8 74 | 336.7 
% del anual 4.4 3.0 10 6.6 102 1 165 | 166 | MS [165 8.1 2.6 PL) US 
El Peaje 

Lluvia media 169 | LS 8.3 (9150 [| SS | 4 [| 23 | Baz [| (607 | 25 9.0 7.0 | 417.2 
% del anual 4.1 2.8 2.0 45 9.4 15.2 | 197 | 140 | 16.7 7.9 2,2 17 | 100.2 
Mexquitic 

Lluvia media 133 | 11.9 6.8 140 | 433 | 722 | 791 | 544 | 70.3 | 33.1 6.9 7.8 | 413.0 
% del anual Si 20) 16 3.4 105 | a | 102 | 2 | 170 8.0 17 19 | 0 
El Terrero 

Lluvia media TUS 3.0 3.0 PMI $89 | Ti | Es | sil [| ¿37 [| 162 8.0 12.6 | 298.6 
% del anual 3.9 10 10 4.1 114 | 158 | 233 | 111 | 16.3 5.4 2.7 4.2 | 100.2 
Xoconoxtle 

Lluvia media 10.8 | 11.9 5.1 22.7 | 379 | 490 | 65.1 | 592 | 73.7 | 28.3 8.3 6.0 | 377.9 
% del anual 2.9 3.1 LS) 6.0 100 | 150 | 72 [ 1837 | 195 ES 222) 16 | 100.0 
V. de Zaragoza 

Lluvia media 8.7 7.8 532 A00 | 833 | $0 [| 47 [| 689 [| 520 [| 25 6.2 4.7 | 308.3 
% del anual 2.8 25 17 6.7 116 | 129 | 158 | 185 | 16.9 7.1 2.0 15 | 100.0 
% anual prom. 3.79 2155 1.33 5.11 | 10.44 | 16.03 | 18.61 | 12.94 | 17.31 | 7.56 2.26 1.92 = 

S del % anual 0.50 | 0.67 | 0.37 | 1.07 | 0.62 | 1.50 | 245 | 144 | 1.22 | 0.91 | 0.31 | 1.09 = 
Cv del % anual | 0.13 | 0.26 | 0.28 | 0.21 | 0.06 | 0.09 | 0.13 | 0.11 | 0.07 | 0.12 | 0.14 | 0.57 = 


se revisó si su lluvia anual mostraba correlación 
con la de la estación Mexquitic, pero no fue así. 
Por lo anterior, se buscaron en el norte del valle 
de San Luis Potosí estaciones pluviométricas 
cercanas, se encontraron dos: Villa Hidalgo y 
Los Pilares; la primera es la más cercana, pero 
fue suspendida en 1986, y la segunda se ubica 
cerca del poblado de Ahualulco. El registro 
disponible de lluvia anual en Los Pilares se 
tiene en la última columna del cuadro 6. 


El diagrama de dispersión entre la lluvia 
anual de las estaciones Los Pilares y El Terrero 
se construyó con veinte parejas de datos, ya 
que en El Terrero faltan ocho años y en Los 
Pilares no se tienen datos en los dos últimos 
años del periodo común. En tal diagrama se 
detectaron cinco años que fueron considerados 
valores dispersos (años 1993, 1995, 1996, 2004 
y 2006); la regresión obtenida con 15 parejas 
de datos tiene un coeficiente de correlación 
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Cuadro 6. Precipitaciones anuales (mm) en el periodo común en las estaciones pluviométricas del valle de 
San Luis Potosí y cercanas. 


AñO Fac. de Los Observ. | Soledad cae El ole Villa de ES 
Ingeniería | Filtros SLP D.G. Terrero Zaragoza | Pilares 
1981 446.7 500.7 444.3 395.9 504.3 492.9 (313.1) 286.3 355.1 474.6 
1982 321.1 326.7 298 195.6 291.4 305.7 (203.4) 197.9 281.3 298.2 
1983 269.3 280.4 317.1 231.0 259.3 285.1 (252.2) 212.9 303.1 376.8 
1984 SES) 361.6 344.6 (307.6) 507.7 564.9 DOYEZ 426.7 (397.4) 772 
1985 477.1 489.1 526.2 347.9 485.3 478.4 388.8 487.7 405.5 437.6 
1986 446.4 456.7 446.2 445.6 403.2 439.7 ANS) 256.9 345.0 421.0 
1987 354.6 364.8 348.9 271.3 322.6 427.0 278.0 143.7 228.5 393.5 
1988 383.6 455.0 438.7 403.8 399.4 576.3 (351.6) 385.7 444.0 536.6 
1989 241.5 239.5 219.0 (216.7) 282.9 270.5 (161.7) 292.9 127.0 231.3 
1990 300: 395.5 341.1 (304.8) PL 420.4 (195.7) 410.0 178.0 285.9 
1991 571.7 528.1 518.2 485.6 476.0 641.2 386.0 709.0 232.9 492.5 
1992 495.0 469.6 414.2 DIZ:0 477.1 480.8 320.0 403.0 181.0 201 
1993 351.6 296.1 318.2 263.9 404.7 (394.6) 459.3 540.0 308.0 382.9 
1994 380.0 484.1 368.7 348.0 (471.1) 386.0 476.0 8960 472.9 
1995 292.4 306.5 278.2 324.7 395.4 386.0 314.0 260.0 337.3 
1996 29 al 295%) 268.0 SIS 304.5 SOÑO 425.0 350.0 AS 
1997 462.6 372.9 151.9 439.0 270.5 (302.8) 466.4 441.5 458.1 
1998 364.3 344.9 (307.5) 8393 160.3 (166.5) 407.0 105.7 238.9 
1999 409.8 306.0 247.3 334.7 163.0 212.0 373.0 184.5 349.9 
2000 1 206.3 276.3 281.9 300.0 125.0 276.0 117.0 219.8 
2001 267.4 312.2 362.8 321.8 499.9 224.3 459.0 145.5 420.9 
2002 289.9 475.8 450.8 556.7 642.2 258.5 418.1 2195 525:2 
2003 470.9 322.5 304.1 319.0 224.7 311.5 332.2 429.0 367.1 
2004 480.9 553.9 537.9 702.7 474.8 245.0 421.0 407.5 607.2 
2005 292.8 213.0 264.3 273.8 358.4 180.0 321.0 355.5 350.9 
2006 510.8 405.7 341.4 556.4 383.0 368.0 436.0 433.1 342.4 
2007 431.3 505.7 399.5 520.2 457.0 343.5 300.5 589.5 578.3 
2008 260.3 323:0 241.3 581.1 590.5 263.1 388.0 376.0 437.5 
2009 497.8 442.7 288.1 533.6 448.4 336.1 288.4 372.0 - 
2010 22 428.3 400.8 563.9 508.4 324.6 484.0 427.5 = 
o Mín 155.2 206.3 151.9 259.3 160.3 125.0 143.7 105.7 219.8 
S Máx 717 ES) 53/19 702.7 642.2 459.3 709.0 589.5 607.2 
E Media 373.8 375.7 329.7 417.2 414.3 288.4 377.9 311.2 400.8 
ss Mediana 360.0 346.9 306.1 401.3 433.3 D99tS 395.5 340.7 388.2 
3 S 98.680 99.934 | 94.575 96.809 114.836 129.778 81.565 113.854 119.199 | 103.653 
5 Cv 0.264 0.253 0.252 0.294 0.275 0.313 0.283 0.301 0.383 0.259 
E Cs -0.026 -0.318 0.092 0.525 0.531 -0.200 -0.107 0.419 -0.025 0.090 
2 Ck 2.625 2.410 2.499 3.005 215 2.720 2.632 4.539 2.764 2.664 
_ Ñ 0.024 0.116 -0.172 0.001 0.012 0.216 0.346 0.164 0.403 -0.109 


Nota: entre paréntesis, los valores estimados por regresión lineal. 
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lineal de 0.713, con un error estándar de la 
estimación de 56.4 mm. Los valores anuales 
estimados en la estación El Terrero se muestran 
entre paréntesis en el cuadro 6; su distribución 
mensual se realizó con base en los promedios 
establecidos en el cuadro 5 y se presentan entre 
paréntesis en el cuadro 3. 

Para completar a la estación Mexquitic, 
también se empleó el registro de lluvia anual de 
Los Pilares, cuyo diagrama de dispersión con 26 
parejas de datos condujo a una regresión lineal, 
cuyo coeficiente de correlación lineal resultó de 
0.659, con un error estándar de la estimación 
de 105.7 mm. Este valor de correlación pudo 
ser mejorado, eliminando datos dispersos de 
arriba y abajo de la recta de regresión, pero ello 
se consideró innecesario. Los valores anuales 
estimados en la estación Mexquitic se muestran 
entre paréntesis en el cuadro 6. Nuevamente, 
su repartición mensual se llevó a cabo con base 
en sus porcentajes, mostrados en el cuadro 5. 

Por último, se dedujo el año faltante en la 
estación Zaragoza con base en el registro anual 
de Xoconoxtle. Su diagrama de dispersión 
mostró dos años dispersos: 1991 y 2007; 
además, nueve puntos se apartan bastante 
de la recta de 45%, indicando valores de lluvia 
anual en Zaragoza muy reducidos en dos 
periodos: 1989-1993 y 1998-2002. Tales datos 
fueron eliminados. Con 18 parejas de datos, la 
regresión mostró un coeficiente de correlación 
lineal de 0.730 y un error estándar de la 
estimación de 45.3 mm. El valor estimado se 
muestra entre paréntesis en el cuadro 6. 


Verificación de requerimientos estadísticos 


Para que los resultados de los análisis pro- 
babilísticos sean válidos teóricamente, la 
serie de datos o muestra utilizada debe cum- 
plir con ciertos requisitos estadísticos: aleato- 
riedad, 
estacionalidad (Campos-Aranda, 2010). Con 
base en los registros de lluvia anual (cuadro 


independencia, homogeneidad y 


6), se probó homogeneidad y estacionalidad 
con dos pruebas básicas: el Test de Helmert 


y el de las Secuencias; además, se aplicaron 
pruebas específicas de persistencia, tendencia, 
cambios en la media y en la varianza (prue- 
bas de Anderson, Sneyers, von Neumann, 
Kendall, Cramer y Bartlett); todas ellas se 
aplicaron con un nivel de confianza del 95%. 
En el cuadro 6, 5, Co, Us , Ek y r, son, res- 
pectivamente, la desviación estándar y los 
coeficientes de variación, asimetría, curtosis y 
de correlación serial de orden 1. Se encontraron 
no homogéneos los registros de la Facultad de 
Ingeniería y de Soledad Diez Gtz., el primero 
según el Test de Helmert y el segundo de 
acuerdo con la prueba de las Secuencias, pero 
ambos por exceso de oscilación. En la estación 
El Terrero se encontró persistencia (r, = 0.346) 
según varias pruebas. Por último, la estación 
Zaragoza resultó no homogénea de acuerdo 
con las pruebas básicas; además, presentó 
persistencia (r, = 0.403), según diversos test. 


Selección por conveniencia estadística de las 
estaciones pluviométricas 


Habiendo quedado por seleccionar una estación 
de las ubicadas en la zona urbana, la que se 
escoge es Los Filtros, pues tal registro resultó 
homogéneo, al igual que el del Observatorio 
SLP; sin embargo, el primero es considerado 
más confiable, ya que ha tenido menos cambio 
de personal, de ubicación y de condiciones 
físicas circundantes. 


Meses de análisis de las duraciones de sequía 


Otra utilidad del cuadro 5 es permitir definir 
los meses correspondientes a los periodos de 
sequía (1, 3, 6 y 12 meses), que serán analizados 
de la manera siguiente. Primero se obtienen 
los promedios mensuales de los porcentajes 
de lluvia, en este caso sólo se promediaron los 
valores correspondientes a las ocho estaciones 
ubicadas dentro del valle, pues se excluyó 
Zaragoza debido a que resultó no homogénea. 
Tales promedios se muestran en negritas. Se 
observa que el mes con más lluvia es julio; 
los tres meses seguidos con mayor lluvia son 
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julio, agosto y septiembre, y los seis meses 
consecutivos con más cantidad de lluvia serán 
de mayo a octubre. 


Obtención de los valores medianos mensuales 


Con base en los registros completos, se pro- 
cesan en cada una de las cinco estaciones 
pluviométricas que serán utilizadas sus re- 
gistros mensuales en el periodo común, para 


obtener los valores medianos tanto muestrales 
como poblaciones, así como otros estadísticos 
centrales (moda y precipitación confiable). 
Los valores citados se obtienen ajustando 
de nueva cuenta la distribución Gamma 
Mixta (Campos-Aranda, 2005) cada mes. Los 
resultados se han concentrado en el cuadro 7, 
mostrando con negritas las magnitudes que 
será utilizadas posteriormente, para obtener 
los déficits de lluvia. 


Cuadro 7. Valores asociados con la tendencia central (mm) en el registro mensual del periodo común 1981-2010 


de las estaciones pluviométricas del Valle de San Luis Potosí que serán procesadas. 


Estación y estadística Ene. | Feb. | Mar. | Abr. | Mayo | Jun. | Jul. | Ago. | Sep. | Oct. | Nov. | Dic. 
Los Filtros 

Núm. de datos nulos 7 8 illo 8 1 0 0 0 0 2 6 ls) 
Mediana muestral 6.5 2.5 0.5 10.8 35.7 | 49.9 70.9 35.5 52.1 25.9 3.8 2.6 
Mediana poblacional 7.8 3.8 0.0 10.9 35.7 47.8 64.1 39.2 52.3 22.6 5102 32 
Moda NE* NE NE NE 23.9 182 | 31.4 16.5 21.0 5.4 0.1 4.5 
Precipitación confiable 0.3 0.0 0.0 0.0 AO) | 247 | ¿82 | 207 | 275 9.4 0.8 0.0 
El Peaje 

Núm. de datos nulos 6 ill 13 7 0 0 0 0 0 2 5 14 
Mediana muestral 7.5 3.0 1.4 10.9 32.8 51.9 69.9 37.1 62.7 | 26.1 4.6 1.5 
Mediana poblacional 8.4 3.3 1.2 11.6 33.1 52.4 68.7 47.2 56.3 PND 5.0 27 
Moda NE NE NE 2.7 19.8 28.4 | 392 | 22.7 | 26.8 4.0 1.4 5.9 
Precipitación confiable 0.7 0.0 0.0 Ml 1613 20 39.6 2D 30.6 9.4 0.0 0.0 
Mexquitic 

Núm. de datos nulos 115 10 16 18) 1 1 0 2 0 1 10 10 
Mediana muestral 3.4 2.8 0.0 4.0 41.6 | 80.5 772 | 493 64.5 28.9 5.0 2.6 
Mediana poblacional 3.2 3.9 0.0 3.6 320 (157-748 59:58 11558 53:08 022357, 4,3 3.2 
Moda NE NE 4.9 1.0 8.8 26.6 14.1 21.3 30.4 2.8 3.3 NE 
Precipitación confiable 0.0 0.0 0.0 0.0 14.9 29.4 28.5 21.0 215 10.0 0.0 0.0 
El Terrero 

Núm. de datos nulos 10 13 145) 4 0 0 0 0 1 3 4 8 
Mediana muestral 6.4 1.0 1.6 8.1 274 | 44.0 594 | 219 | 446 12.3 4.9 6.8 
Mediana poblacional 4.9 0.8 1.3 8.5 20 [| sE | Eno | 222 | “iso | 120 5.6 5.8 
Moda 0.1 0.1 1.9 0.9 16.9 6.2 42.1 5.0 31.0 4.7 1.2 NE 
Precipitación confiable 0.0 0.0 0.0 zo) 16.9 RN IT) ESA 2610) Srl 1.8 0.0 
Xoconoxtle 

Núm. de datos nulos 9 12 12 6 1 0 0 0 1 S 7 16 
Mediana muestral 2.0 11 1.0 10.7 | 345 | 415 58.5 | 430 | 640 | 23.5 4.5 0.0 
Mediana poblacional 4.0 14 1.4 10.7 30.1 SO 53.0 40.8 60.7 21.0 4.9 0.0 
Moda NE NE NE NE 13.5 3.4 26.5 NE 33.5 10.2 0.5 1.8 
Precipitación confiable 0.0 0.0 0.0 0.7 115522 153 | 222 TAO 8223: 6.9 0.4 0.0 


* NE: no existe. 
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Interpolación espacial de registros 


Más del 90% del área del valle de San Luis 
Potosí queda comprendida en la carta Hi- 
drológica de Aguas Superficiales F14-4, 
denominada San Luis Potosí. La información 
de esta carta se utilizó para trazar el 
parteaguas de tal cuenca cerrada en la carta 
topográfica (INEGI, 1996) de igual clave, 
ambas impresas y con escala de 1:250 000. 
El área del valle asciende a los 1 858 km?. En 
la figura 3 se observa que la forma del valle 
se puede aproximar a la de un rectángulo, 
cuyos lados se dividieron en cinco tramos, 
para formar una malla de 25 rectángulos, que 
son casi cuadrados y por ello, en adelante así 
se llamarán. Cada rectángulo tiene una altura 
aproximada de 8.5 km y una base de 9.25 km, 
por lo cual su área es de 78.625 km?. Con 
tales dimensiones, el valle de San Luis Potosí 
quedó representado por 23 cuadros (97.3% 
del área), como se muestra en la figura 3. 

En seguida se localizó el centro de 
gravedad de cada rectángulo, por el punto de 
cruce de sus bisectrices. Después se midieron 
manualmente en la carta topográfica F14-4 
impresa las distancias (d,) en milímetros, entre 
cada centro de un cuadrado y cada una de 
las cinco estaciones pluviométricas que serán 


utilizadas. Los resultados se han concentrado 
en el cuadro 8. Con base en las distancias d, se 
calculan los cinco factores de ponderación para 
cada cuadrado (a) y con ellos se ponderan 
los registros mensuales, para estimar la serie 
hipotética. Lo mismo se realiza con los valores 
medianos del cuadro 7, a fin de estimar las 
magnitudes mensuales hipotéticas de cada 
cuadro, cuyos valores anuales se indican en la 
última columna del cuadro 8. 


Cálculo de los déficits anuales de lluvia 


Estimadas en cada cuadro, la serie hipotética 
de lluvia mensual en el periodo común y 
los valores ponderados de las medianas 
mensuales, se calcularon los treinta déficits 
anuales de la lluvia para duraciones de 12, 6, 
3 y 1 meses. Tales magnitudes son la suma de 
sólo los valores negativos obtenidos en cada 
año del periodo común, en todo el año, en el 
lapso de mayo a octubre, durante julio, agosto 
y septiembre, y en julio. Por limitaciones de 
espacio, en adelante sólo se expondrán los 
resultados para las duraciones de sequía 
de 12 y 6 meses. Calculados de la manera 
descrita, los déficits son negativos, pero se 
indican y procesan probabilísticamente sin 
el signo menos. Mayores detalles sobre este 


Cuadro 8. Distancias (mm) entre cada centro de gravedad del cuadrado de la malla y la estación pluviométrica indicada. 


PMA es la precipitación media anual obtenida en el periodo común y MED la mediana anual estimada por ponderación. 


Núm. del Estación número: PMA | MED | Núm. del Estación número: PMA | MED 
cuadro 1 2 3 4 5 (mm) | (mm) | cuadro 1 2 3 4 5 (mm) | (mm) 
il 7 | 0 | E 29 (170 | 632 | 205 13 is [| 169 [187 [| ES 85 | 364.7 | 267.9 
2 117 | 156 | 97 13 | 157 | 321.6 | 243.8 14 37 10 78 | 158 | 157 | 405.9 | 302.2 
$ 144 | 186 | 133 | 25 | 152 | 333.0 | 250.4 il; 5 46 TO STA 23 SO 2.59 29055 = 
4 65 66 19 | 127 | 191 | 399.3 | 289.0 16 41 83 9% | 10 91 | 3845 | 282.4 2 
5 55 80 21 91 | 163 | 394.0 | 285.9 il 77 | 0 || 9 64 | 376.8 | 275.1 3 
6 67 | 106 | 57 58 | 140 | 373.5 | 274.8 18 113 | 157 | 151 | 101 51 | 373.1 | 271.0 7 
7 9 | 137 | 93 29 [| 165 | 3502 | 2205 19 58 sy | M2 | 11 | 150 | $056 | 2925 El 
8 126 | 170 | 131 39 | 118 | 348.1 | 259.0 20 46 59 | 112 | 158 | 113 | 390.3 | 287.3 E 
9 41 88 46 | 139 | 172 | 397.3 | 292.4 21 61 91 | 124 | 141 77 | 384.4 | 281.3 E 
10 21 56 47 | 108 | 141 | 391.7 | 288.2 22 89 | 125 | 144 | 133 | 42 | 379.3 | 274.6 E 
11 45 89 71 81 | 113 | 381.0 | 280.1 23 1229 LT 11709 15 OA ES7715 [4268:2 S 
12 79 | 124 | 112 66 | 93 | 368.8 | 271.5 Prom. - - - - - 376.0 | 276.4 . 
E 
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procedimiento se pueden obtener en Campos- 
Aranda (2012). 


Procesamiento probabilístico de los déficits 
anuales 


Cada una de las cuatro series de treinta valores 
del déficit anual, obtenidas en cada uno de los 
23 cuadrados de la malla (ver figura 3), fueron 
procesadas  probabilísticamente, ajustando 
la distribución General de Valores Extremos 
(GVE), a través de los métodos de sextiles y 
de momentos L (Campos-Aranda, 2010). Se 


emplearon sólo estos dos métodos de ajuste, 


pues en algunas de las series citadas ocurrieron 
ceros. En el cuadro 9 se han concentrado las 
predicciones obtenidas para los déficits anuales 
de duración de 12 meses, en los seis periodos 
de retorno que fueron establecidos para las 
curvas SAF. 

En cambio, en el cuadro 10 se muestran las 
predicciones de los déficits anuales de duración 
de seis meses, sólo en los periodos de retorno de 
5, 10, 25 y 50 años. Las predicciones indicadas 
en los cuadros 9 y 10 proceden del método 
de estimación de parámetros de la GVE, que 
condujo al menor error estándar de ajuste 
(Kite, 1977), por lo general, el de momentos L. 


Cuadro 9. Predicciones (mm) de los déficits anuales (duración = 12 meses) estimados en cada centro de gravedad 


de un cuadro de la malla (Local) y sus valores máximos por áreas contiguas (VMAC), en el valle de San Luis Potosí. 


Periodo de retorno en años 
Núm. del | EEA* 
A A 5 10 25 50 75 100 

Local | VMAC| Local [VMAC| Local [|VMAC| Local [VMAC| Local [VMAC| Local [| VMAC 

1 5.2 70.8 92.3 91.6 | 113.0 | 120.1 | 143.4 | 142.9 | 169.1 | 156.8 | 184.8 | 167.1 | 196.3 

2 8.1 65.7 90.2 82.6 | 112.2 | 103.8 | 141.1 | 119.3 | 167.1 | 128.3 | 183.0 | 134.6 | 194.8 

3 7.2 64.3 88.3 82.2 | 110.5 | 106.4 | 140.1 | 125.6 | 164.9 | 137.3 | 180.2 | 145.8 | 191.4 

4 5.9 89.8 88.7 113.0 | 111.1 | 141.3 | 139.4 | 161.7 | 164.1 | 173.2 | 178.4 | 181.3 | 189.1 

5 57 85.5 88.0 108.5 | 110.6 | 137.1 | 139.7 | 158.1 187.6 

6 ll 73.4 87.5 94.8 | 109.9 | 124.2 | 139.3 | 148.0 186.3 

7 6.1 65.9 86.8 85.0 | 109.3 | 1119 | 138.5 | 134.1 185.3 

8 6.1 65.0 86.1 a | tas | 10s | es | 1808 184.3 

9 5.2 84.6 85.2 | 107.1 | 107.5 | 137.1 | 136.9 | 160.5 183.7 

10 6.3 83.1 844 | 106.4 | 106.6 | 138.8 | 136.1 | 165.0 182.7 

11 5.4 76.4 83.8 98.2 | 105.8 | 128.5 | 135.2 | 153.0 181.8 

12 5.4 71.0 84.0 90.8 | 105.7 | 118.1 | 1343 | 140.2 181.1 

13 5.2 69.2 83.4 88.2 | 105.1 | 1144 | 133.3 | 135.5 180.0 

14 5.8 928 82.7 | 113.0 | 104.4 | 138.1 | 133.0 | 156.1 179.0 

a 15 6.7 88.0 82.2 111.4 | 103.7 | 143.4 | 132.2 | 169.1 177.9 
2 16 27) 78.2 81.7 99.8 | 103.1 | 129.8 | 131.4 | 154.0 176.7 
3 17 4.8 74.1 81.0 94.1 | 1024 | 121.2 | 130.7 | 142.7 175.8 
a 18 4.8 74.4 80.5 Sa Wi | US | 1299 | 1661 174.8 
E 19 6.1 84.5 79.9 105.9 | 101.0 | 133.9 | 129.1 | 155.5 174.0 
E 20 6.2 80.7 2 | 125 | 1002 | 1625 [| 1ES | 1887 173.0 
hal 21 5.3 77.4 78.5 98.3 99.4 | 126.6 | 127.4 | 149.3 1721 
3 20) Sul 78.2 TS) 978 98.7 | 121.7 | 126.4 | 140.2 170.9 
a 23 6.5 86.0 77.3 105.3 | 97.9 128.7 | 125.4 | 145.3 | 146.9 | 154.7 | 159.9 | 161.2 | 169.3 


* Error estándar de ajuste. 
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Cuadro 10. Predicciones (mm) de los déficits de seis meses de duración estimados en cada centro de gravedad de 
un cuadro de la malla (Local) y sus valores máximos por áreas contiguas (VMAC), en el valle de San Luis Potosí. 


Núm. del EEA* Periodo de retorno en años 
cuadro (mm) Sl E a a 
Local VMAC Local VMAC Local VMAC Local VMAC 
1 5.0 62.5 78.9 81.8 100.7 107.9 128.6 128.7 153.3 
2 6.6 56.2 IDA PES 98.6 93.6 1259 109.9 ISI4S 
3 6.2 54.8 76.0 71.8 97.3 95.5 125.4 114.8 149.1 
4 57 78.4 76.6 100.7 OEA 128.6 126.1 149.1 149.1 
5] 6.5 75.5 76.4 96.5 97.6 123.2 125.6 143.1 147.9 
6 50) 64.8 75.8 84.6 97.0 11.6 1511 SS 147.1 
dl 5.0 57.2 75.2 75.3 96.5 101.0 124.4 122.3 146.3 
8 Sl 56.4 74.6 74.0 95.8 99.1 123.8 MES 145.7 
9 5.6 73.6 73.8 94.6 94.9 122.5 123.1 144.2 145.1 
10 6.3 71.8 73.1 93.7 94.1 124.4 para 149.7 144.5 
1 5.4 66.5 72.6 87.1 93.5 115.7 121.6 139.1 143.9 
12 4.9 62.0 Ya 80.9 5) 107.1 120.9 128.4 143.0 
13 4.9 60.5 72.6 78.8 93.4 104.4 121.1 125.4 142.2 
14 5.9 78.9 22! 98.9 27 187 120.4 141.8 142.2 
15 LZ 75.4 71.6 97.5 92.2 128.2 119.7 153.3 141.5 
16 50) 67.6 ZN EA 88.2 91.6 1407/22 119.1 141.1 140.7 
17 5.1 64.7 70.7 83.9 91.0 110.4 118.4 131.8 140.0 
18 4.8 64.9 70.2 83.5 90.5 108.8 117.8 129.0 139.4 
19 6.1 73.1 69.7 93.4 89.9 120.3 117.1 141.1 138.6 
20 6.3 69.9 69.1 90.7 89.1 ENS 116.3 142.9 137.8 
21 5.8 67.1 68.4 87.1 88.4 114.7 115.4 137.2 137.0 
20 2 67.8 67.9 86.8 87.7 112.4 114.5 132.6 136.0 
23 6.8 78.8 67.3 98.6 87.0 123.7 113.6 142.5 134.8 


* Error estándar de ajuste. 
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Figura 4. Curvas severidad-área-frecuencia de duración de 12 meses del valle de San Luis Potosí. 
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Figura 5. Curvas severidad-área-frecuencia de duración de seis meses del valle de San Luis Potosí. 


Integración de las curvas SAF 


Las predicciones de cada uno de los periodos 
de retorno procesados (cuadro 9) se vaciaron 
en su cuadro respectivo y se procedió a obtener 
en tales mapas, uno para cada intervalo 
de recurrencia, los valores máximos por 
áreas contiguas, los cuales también se han 
concentrado en el cuadro 9 para la sequía 
con duración de 12 meses. Con tales valores 
se construyen las curvas SAE las cuales se 
muestran en la figura 4. 

Siguiendo un proceso idéntico, se formaron 
los mapas de predicciones, según datos del 
cuadro 10, para obtener los valores máximos por 
áreas contiguas, para la sequía de duración de 
seis meses; con estas magnitudes se elaboraron 
las curvas SAF respectivas, mismas que se 
presentan en la figura 5. Mayores detalles de tal 
procedimiento se pueden consultar en Mishra 
y Singh, 2008. 


Análisis de las sequías históricas 
Durante las etapas de estimación de los déficits 


anuales y su análisis probabilístico, se detectó 
que los años de 1982, 1989, 1995, 2000 y 2005, 


presentaban los déficits más grandes y por 
ello son años de sequía. Para estos cinco años 
se estimaron por ponderación en cada cuadro 
de la malla sus registros hipotéticos mensuales, 
así como sus valores medianos mensuales. 
Formados tales registros, se obtuvieron, en 
cada cuadro de la malla, los déficits anuales 
de duración de 12, 6, 3 y 1 meses. Por ejemplo, 
en el cuadro 11 se exponen los relativos a la 
duración de 12 meses en los cinco años de 
sequía definidos previamente. También en 
este cuadro se muestran los déficits anuales de 
duración de seis meses, pero sólo en los años 
de 1982 y 2000. 

De nuevo se construyen mapas de cada año 
de sequía, con los valores de los déficits anuales 
por cuadro y se obtienen los valores máximos 
por áreas contiguas, los cuales se dibujan en 
las curvas SAF de la duración correspondiente, 
para observar la frecuencia o periodo de 
retorno que le corresponde a tal sequía, así 
como su evolución espacial dentro del valle 
de San Luis Potosí. Las sequías del cuadro 11 
se tienen dibujadas con línea punteada en las 
figuras 4 y 5. 

En la figura 4 se observa que todos los 
años de sequía siguen de manera bastante 
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Cuadro 11. Déficits (mm) calculados en cada centro de gravedad de un cuadro de la malla (Local) y sus valores 


máximos por áreas contiguas (MAC), en las siete sequías indicadas del valle de San Luis Potosí. 


Sequías analizadas de duración de 12 meses Sequías de duración de seis meses 
Núm. del 
e 1982 1989 1995 2000 2005 1982 2000 
Local [| VMAC| Local | VMAC | Local | VMAC| Local | VMAC| Local | VMAC| Local | VMAC| Local | VMAC 
1 23.0 113.6 109.8 138.9 93.8 140.5 144.9 154.2 88.0 99.6 JT 104.6 131.6 139.8 
2 43.1 109.3 74.2 137.1 80.3 137.8 142.7 154.0 83.9 99.0 43.1 101.0 128.5 1897 
3 52.0 107.5 83.2 136.3 86.5 134.8 141.0 152.7 82.0 98.0 51,6 99.4 127.7 137.7 
4 101.1 106.5 138.5 135,9 98.0 132.5 154.2 152.3 99.6 98.3 9382 98.6 139.8 136.5 
S 95.6 105.4 134.6 136.4 97.1 130.4 153.7 150.9 98.3 96.4 90.1 97.9 139.6 135.5 
6 81.2 104.6 116.8 136.0 100.3 128.7 145.1 149.9 88.1 95.0 76,9 371 131.6 134.9 
Ed 64.3 103.9 97,2 135.3 94.0 127.4 141.4 149.1 82.9 94.0 62.2 96.5 128.4 134.2 
8 64.5 102.9 95.3 134.5 94.6 126.5 138.8 148.3 79.6 93.1 62.6 997 126.7 193.5 
9 103.4 101.9 133.8 133.5 119.6 125.4 141.7 147.5 88.9 92.2 96.2 94.9 125.1 132.9 
10 100.6 101.1 135.2 132.6 121.8 124.8 150.2 147.0 96.2 91.4 93.8 94.2 133.6 132.4 
11 89.9 100.8 123.1 131.8 112.6 123.8 144.2 146.4 86.8 90.7 84.4 93.9 130.1 132.0 
12 79.3 101.5 112.3 130.8 107.2 122.8 138.2 146.0 78.9 90.0 dl 94.7 126.3 181,5 
13 79.6 100.7 108.6 129.8 104.6 122.0 134.9 145.5 74.4 89.4 75.9 94.1 124.6 131.1 
14 113.6 999 134.9 129.0 140.5 120.9 126.5 145.0 79.3 88.9 104.6 93.4 106.8 130.7 
15 104.9 99.2 138.9 128.2 135.1 119.8 151.1 144.5 99.0 88.2 97.3 92.7 132.9 130.3 
16 94.3 98.0 125.8 127.3 120.2 118.5 141.7 144.0 84.6 87.7 88.2 91,7 127.8 1299 
17 88.7 96.9 118.6 126.5 113.5 117.4 133.6 143.5 72.6 87.2 83.5 90.8 123.7 129.5 
18 ST 96.0 sd 125.5 111.8 116.3 127,5 143.0 64.2 86.6 86.8 89.9 120.6 1292 
19 104.0 94.9 130.9 124.6 128.8 115.1 134.8 142.5 82.0 86.0 96.4 89.0 118.2 128.8 
20 207 93.4 129.0 133 125.6 114.1 1592 141.7 84.1 85.3 92.8 87.7 124.1 128.3 
21 94.0 92.0 124.6 121.9 120.1 113.1 135.8 141.0 76.9 84.3 87.8 86.4 123,9 127.8 
2 97 90.2 120.0 120.2 117.0 111.9 125.8 140.3 62.0 83.3 91.4 84.8 119.4 127.1 
23 109.1 88.1 117.6 118.2 119,3 110.5 114.2 139,2 44.3 81.6 103.3 83.1 111.9 126.2 


aproximada el comportamiento obtenido en 
las curvas SAFE lo cual significa que en el valle 
de San Luis Potosí, las sequías presentan el 
mismo patrón de evolución espacial. De los 
cinco años con sequías históricas, la del 2005 
fue la más leve, cuyo periodo de retorno es del 
orden de seis años; en cambio, la sequía del 
año 2000 fue la más severa, con un intervalo 
de recurrencia del orden de los cuarenta años. 
Por otro lado, al estudiar la primera y última 
sequías (1982 y 2000) del periodo común 
analizado, en las curvas SAF de duración de 
seis meses de la figura 5, también se observa 
que siguen el comportamiento general de las 
curvas y que ahora la sequía de 1982 tiene un 
periodo de retorno de diez años. 


Modificaciones al procedimiento 
aplicado 


Existe un refinamiento que puede ser intro- 
ducido para mejorar la interpolación de los 
registros en los centros de gravedad de los 
cuadrados de la malla y consiste en utilizar la 
técnica del Kriging (Metcalfe, 1997). También 
se visualiza que todo el procedimiento descrito 
puede ser computarizado o bien desarrollado 
en un sistema de información geográfica, 
para su automatización. Lógicamente, para la 
cuantificación de las sequías se pueden utilizar 
otros índices, como el de la precipitación 
estandarizada (Velasco-Velasco, 2002; Mishra 
y Singh, 2008; González-Camacho et al., 2011). 
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Por otra parte, el valle de San Luis Potosí 
es en realidad una cuenca pequeña dentro 
del análisis espacial de sequías, y por ello, los 
mapas construidos con las predicciones de los 
déficits anuales resultaron sumamente simples, 
mostrando los valores máximos en la zona 
montañosa del oeste del valle y los mínimos 
al norte. Sin embargo, en cuencas mayores, es 
muy probable que existan algunos puntos o 
cuadros con valores máximos y entonces habrá 
que estudiar diversas “evoluciones espaciales” 
de la sequía para obtener sus valores máximos 
contiguos, de manera similar a como se realizan 
los análisis de precipitación-área-duración, 
al procesar los variados “ojos de tormenta” 
encontrados en las curvas isoyetas (Pérez- 
Méndez y Campos-Aranda, 2000). 


Conclusiones 


Las curvas de severidad-área-frecuencia (SAF) 
se pueden construir en una cuenca o región 
previamente definida, para diversas duraciones 
de la sequía, por ejemplo 12, 6, 3 y 1 meses. 
Tales curvas mostrarán la intensidad de las 
sequías de acuerdo con el índice seleccionado 
y utilizado en su estimación; además, carac- 
terizarán a las sequías en relación con el 
modelo de comportamiento espacial establecido 
por las predicciones procedentes del índice 
empleado. Las curvas SAF de una cuenca o 
región permiten su comparación con otras 
para apreciar sus diferencias, sean éstas en 
intensidad o en comportamiento espacial. 
Estimadas las curvas SAFE, permitirán ca- 
racterizar los años con sequías en tal cuenca 
o región, para observar cómo cambia su seve- 
ridad a lo largo del área y dentro del periodo 
de retorno (frecuencia). Lo anterior ocurre 
al dibujar en las curvas SAF respectivas, los 
valores máximos contiguos del año con sequía 
histórica y observar que éstos las “cruzan hacia 
abajo”, indicando con ello que tal año comenzó 
con sequías muy intensas en áreas reducidas 
y después disminuyó en severidad, al abarcar 
más extensión, o por el contrario, empezó 
con baja intensidad y conforme fue creciendo 


en área, la sequía fue cada vez más severa. 
Esto último no fue observado en las sequías 
históricas del valle de San Luis Potosí, quizás 
debido al tamaño reducido de su cuenca. 

El procedimiento sugerido para estimar las 
curvas SAF de las sequías de una cuenca o re- 
gión geográfica es bastante sencillo, y suscep- 
tible de mejorarse y automatizarse compu- 
tacionalmente. Aunque resulta tedioso en la 
parte de estimación de los valores anuales y 
mensuales faltantes, la obtención y definición 
de un amplio periodo común en las estaciones 
pluviométricas utilizadas conducirá a 
predicciones más confiables y, por lo tanto, 
tales curvas SAF tendrán una mayor garantía 
o exactitud. Lógicamente, el uso de una red 
amplia y uniforme de pluviómetros conduce 
a unas curvas SAF confiables. Por lo anterior, 
estos dos aspectos del procedimiento requieren 
atención especial. 
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Abstract 


CAMPOS-ARANDA, D.E Meteorological drought characterization using severity-area- 
frequency curves, in San Luis Potosi Valley of Mexico. Water Technology and Sciences (in 
Spanish). Vol. IV, No. 3, July-August, 2013, pp. 165-186. 


This work is composed of three sections. The first briefly discusses general ideas about 
droughts, including definitions and types. The second section provides a detailed presentation 
of the procedure suggested to estimate severity-area-frequency curves (SAF). The third part 
describes a numerical application for the valley of the city of San Luis Potosi , Mexico. The 
proposed sequence consists of five basic stages: (1) compilation and selection of records to 
be used, (2) spatial interpolation of records, (3) calculation of annual rainfall deficits, (4) 
probabilistic processing of annual deficits and (5) construction of SAF curves. The temporal 
and spatial characteristics of the droughts are essential to understand their behavior and 
evolution in a watershed or specific region. The procedure for a regional estimation of the 
SAF curves was established using one of the simplest indices to characterize droughts and a 
basic spatial interpolation technique. This included the use of the generalized extreme values 
distribution to obtain predictions that define the frequency of these curves. Lastly, conclusions 
are formulated which point out the usefulness of SAF curves to compare droughts between 
regions and to interpret years of historical drought. The simplicity of the proposed method is 
also shown and the first stage is indicated as the most critical. 


Keywords: monthly precipitation, median values, homogeneity, spatial interpolation, 
predictions, GEV distribution. 
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